Guidance and Control in a Josephson Charge Qubit 
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In this paper we propose a control strategy based on a classical guidance law and consider its use 
for an example system: a Josephson charge qubit. We demonstrate how the guidance law can be 
used to attain a desired qubit state using the standard qubit control fields. 
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INTRODUCTION 



The practical operation of devices for quantum informa- 
tion processing is dependent upon the ability to con- 
trol the behaviour of the component qubits via exter- 
nal classical control fields. As in classical devices, the 
control/bias fields allow the operator to define the dy- 
namical characteristics of the system. Bias fields contain 
noise that will couple into the device and can ultimately 
limit the coherent evolution of a quantum information 
processing system. In classical systems, a feedback con- 
trol loop is often used to reduce the effect of such envi- 
ronmental noise or other unforseen perturbations on the 
evolution of system. The control loop compares the de- 
sired behaviour with the actual behaviour of the system 
and aims to minimise the error between the two. The 
problem for quantum control is that a feedback loop re- 
quires some form of measurement to be made, and this 
measurement will often adversely affect the coherence of 
the quantum evolution. Several groups have suggested 
methods to overcome this problem, using techniques de- 
veloped (mainly in quantum optics) to describe 'weak' 
measurements. These measurements can be used to ob- 
tain information about a quantum system over a period 
of time whilst minimising the adverse effects of the mea- 
surement interaction (see reference Q for a recent re- 
view and a description of the relationship between clas- 
sical control and quantum control). Closed-loop tech- 
niques fall into two main categories: Markovian feedback 
0,0,0 and Bayesian (or optimal) feedback 0,0,0. The 
first method uses the results of measurements to directly 
alter the external control fields applied to the system. 
The second builds an estimate of the system state over a 
number of measurements. Although they were developed 
in quantum optics, these techniques and related analysis 
have recently been applied to the control of solid-state 
qubits 00. 

This paper deals with an associated problem, that of 
generalising the techniques of classical guidance (see for 
example [2j) to the operation and manipulation of qubits. 
The main conceptual difference between guidance and 
control is one of timeliness. In control systems, the de- 



sired state of the system (classical or quantum) may be 
static or change with time, but there is always an error 
between the actual state of the system and the desired 
state. A control is applied to remove this error signal. 
In guidance systems, the evolution of the system is not 
as important as the final state. The controls are ap- 
plied throughout the evolution to ensure that the system 
reaches the desired state at the desired time. In 01 
Bouten et al. have addressed this problem implicitly, by 
using dynamical programming to solve an optimal control 
problem by minimising the controls (which define a 'cost' 
function) applied over the time available (t max ). Experi- 
ence with classical guidance techniques shows that, whilst 
such algorithms may give a minimum cost solution, op- 
timal control guidance can be difficult to implement and 
simpler guidance laws often provide sufficient accuracy 
with significantly simpler guidance-control systems 
The most commonly used guidance law is referred to as 
proportional navigation, which is used in a wide variety of 
aerospace guidance systems (autopilots, guided missiles, 
etc.). Several variants of proportional navigation exist, 
but - in its the most general form - it can be written as 

0, 



a r = N 



, (ZEM) 



' go 



where a c is the control (acceleration command) that 
should be applied to the system, N' is a constant (called 
the 'navigation constant') which determines the strength 
of the commands, t go is the time to go until the objec- 
tive (t max > t go > 0), and ZEM is the 'zero effort miss' 
(that is, the distance between the desired state - the in- 
tercept point - and the predicted state if no more con- 
trols are applied). In classical guidance, an intercept is 
assured as long as N' > 2 and the accelerations com- 
manded are achievable. In practice, N' is normally in 
the range AT' = 4 — > 6, so that the controls immediately 
prior to intercept are minimised. 

Proportional navigation guidance is not optimal in the 
sense of minimising the controls applied, but is generally 
easier to implement in a practical control system and the 
controls that need to be applied to the system tend to 
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have a lower bandwidth than those generated by more 
sophisticated algorithms. In this paper, we will gener- 
alise this classical guidance law to the problem of control- 
ling a solid-state qubit. We will show some examples of 
the behaviour predicted for an (open-loop) proportional 
guidance law applied to a Josephson charge qubit (e.g. 
[Ill Il2j p. (Open-loop control, i.e. without feedback, has 
been studied in the context of atomic and quantu m op - 
tics but not in the same form as that presented here 
In particular, we consider the robustness of the guidance 
law to noise in the bias fields and the introduction of 
a first-order time delay (low pass filter) into the control 
system, to investigate the effect of restricting the band- 
width of the control signal. Consideration is also given 
to the introduction of a simple measurement interaction 
and feedback control loop. 



JOSEPHSON CHARGE QUBIT 

The qubit studied in this paper is a commonly used 
(idealised) model for a standard experimental configu- 
ration. It consists of a superconducting island (a Cooper 
pair box) coupled to an external circuit via two paral- 
lel Josephson junctions 

[HE3 (see Figure 1). The qubit 
has two main control fields, a voltage bias (V x ) to control 
the energy of the charge states, and a magnetic flux (Q x ) 
to control the tunnelling of electron pairs between the 
box and the external circuit. The two parallel Josephson 
junctions form a current loop and applying a magnetic 
field through this loop acts so as to modulate the tun- 
nelling rate onto and off the Cooper pair box. (In this 
paper, we assume that the two Josephson junctions are 
identical for simplicity. In practice, there will be small 
variations in the tunnelling rates for each junction in any 
experimental system and it might be necessary to char- 
acterise these differences in a real system). The effect 
of modulating the effective tunnelling frequency on the 
qubit energy levels is shown in Figl(b). We shall use cir- 
cuit parameter values based on the experimental values 
given in 0, E3 to ensure that the circuit parameters 
are realiseable. In most experiments that have been re- 
ported using such systems (as well as in other supercon- 
ducting qubit experiments based on persistent current 
devices [U4), excited states are generated in the qubit by 
applying an additional field, a time-dependent microwave 
drive field. In this paper, we do not use an additional (ex- 
ternal) microwave drive which reduces the complexity of 
the control system. This point it discussed in more detail 
below. 

The Hamiltonian for the (two-state) qubit can be writ- 
ten in the charge basis representation as [HI fl^ , 
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FIG. 1: (Color online) (a) Schematic diagram of qubit circuit, 
showing capacitances Ci + C2 = C; (b) Energy levels for 
qubit as a function of V x , using parameters given in the text 
and showing energies for nominal bias point <& X(I — 0.25$o 
(solid lines) and for extremes of the control modulation field 
$ xo = 0.25 ± 0.05*o (dashed lines). 



where the basis states are zero excess pairs (|0)) and one 
excess pair (|1)) on the Cooper pair box, C is the net 
capacitance of the island/box (in this case, we take C — 
6 x 10~ 16 F), and v is the tunnelling (angular) frequency 
of the Josephson junction. The maximum Josephson tun- 
nelling frequency is taken to be v/2-k ~ 12.9 GHz, al- 
though the effective Josephson frequency at the nominal 
bias point & x = <I>o/4 is approximately 9.1 GHz (see be- 
low), in line with the parameters given in where 
$0 = h/2e = 2 x 10~ 15 Wb is the superconducting flux 
quantum. For the purposes of this paper, we take the 
charge basis to be the computational basis for the qubit. 
The energy eigenstates are functions of the bias fields. 
Although we we will consider transitions between energy 
eigenstates under the action of the guidance/control, this 
is not necessary. As discussed below, we use the ground 
state as the initial state for convenience, since it is as- 
sumed that the system will relax to this state after some 
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period of time, under the action of whatever dissipation 
processes are present in the system. 

This Hamiltonian (or any other 2x2 Hermitian ma- 
trix) may be decomposed into four components, which 
correspond to a constant multiplied by the identity ma- 
trix (J) or one of the three Pauli matrices 



1 

1 



-i 

1 



We write the bias voltage as V x — e/C - 
decompose the Hamiltonian as, 



1 

-1 

AV X , and then 



H 



Hjl + H x a x + HyUy 

( C(AV x f e 2 
\ 2 2C 



I 



H z a z 

- e(AV x )a z 



hv 



(2) 



From this decomposition, it is easy see that there is no 
<j y term in the basic Hamiltonian. This would be the 
term that would normally be responsible for exciting the 
qubit into the excited state and has been the most com- 
mon control coupling to be studied in quantum optics 
0, IE IHl El- To generate a a y term, it is necessary 
to apply a time-dependent field, since a x and a z do not 
commute: [<7 zi a x ] = 2ia y . Normally, in quantum optics, 
a laser is used to pump a qubit into an excited state, or in 
solid-state experiments (such as those described in refer- 
ences E3 an d EH) an external microwave source is used. 
In an experimental system, the underlying Hamiltonian 
may not be exactly what is predicted by the idealised 
model used here, but a number of techniques have been 
proposed to allow the deviations to be characterised E3 • 
The use of an external microwave source in solid-state 
is not ideal for large scale systems because of potential 
problems in isolating qubits from drives applied to neigh- 
bouring devices. Because of this, we restrict ourselves to 
controls that arise from time-dependent bias fields AV X 
and $> x and consider the effect of limiting the bandwidth 
of these fields in a later section. 

The general representation for a qubit state is, 



COS|- 



sm|- 



11) 



(3) 



where 9 £ [0, it] and <fr € [0, 2ir], which can be written as 
a (pure state) density matrix p = \tp)(tj)\, 



cos 2 (I) 
.(|)sin(f)e^ 



(§)sm(j)e-* 
sin 2 (f) 



(4) 



However, the most convenient representation for the pur- 
poses of this paper is the Bloch sphere representation |17| , 
where the two angles 9 and <p represent angles on a unit 
sphere, defined in a three-dimensional space by, 



(5) 




The different components (a x , a y and a z ) present in the 
Hamiltonian represent rotations in this three-dimensional 
Bloch space (about the X, Y and Z axes respectively). 
The fact that there is no a y term in the basic Hamilto- 
nian is not a problem, because it is possible to reach any 
point on the Bloch sphere from any other by successive 
rotations about any two (non-parallel) axes. The guid- 
ance algorithm simply governs the size of the rotations 
that are to be applied to achieve the objective. 



PROPORTIONAL GUIDANCE 

The classical proportional navigation algorithm predicts 
the expected miss distance if no control is applied (the 
'Zero Effort Miss' or ZEM), a quantum analogue for 
the qubit can be developed in the similar manner. The 
evolution in the absence of controls is described by the 
basic Hamiltonian given in equation (2) and the time 
evolution of the wavefunction and (pure state) density 
matrix is governed by a unitary evolution operator, 



U{t) = exp 



iHnt 



The main difference between proportional navigation on 
a sphere and proportional navigation in three Euclidean 
dimensions is that the rotations generated by this unitary 
matrix and the rotations required to move the estimated 
final state onto the desired state do not commute. The 
rotation required at the end point will not produce the 
same effect if it is applied earlier in the trajectory. Be- 
cause of this, we need to retrodict where the desired state 
should have been at the earlier time, if it is to end up 
at the desired state under the free evolution given by 
Hq. The ZEM in this case is the two angles (9zem and 
4>zem) which separate the current estimated state (on 
the Bloch sphere) from the point where the desired state 
would have to be at the current time. So, we calculate 



Pd{tgo) 



U\t go ) Pd U{t go ) 
iHptgo 



exp 



Pd exp 



iH t qo 



(6) 



where pd is the desired final state (or 'target state'). 
From this density matrix, we can calculate the two angles 
9d(t go ) and <pd(t go ) which define the retrodicted state. 
The ZEM angles are then given by, 



^ZEM — 9d(t go ) — 9 
^ZEM = 4>d{tgo) - 4> 



(7) 



where 9 and <p represent the current state (allowing for 
the periodicity of the angles). The controls that need to 
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be applied are 



M 
dt 
(hp, 
~~dt 



^£ = N '- 



>ZEM) 



zem) sin ( 



tgo 



(8) 



where the controls are angular velocities rather than ac- 
celerations because the Bloch equations are first-order 
differential equations rather than second-order clas- 
sical dynamics, and the sin (9 term arises because the dif- 
ferences in 4>zem near the poles of the Bloch sphere need 
to account for the curvature of the sphere. 

The controls can be equated to an equivalent Hamil- 
tonian by integrating over a small time interval (St) and 
using the fact that a cr x -Hamiltonian generates rotations 
about the X-axis and a cr z -Hamiltonian generates rota- 
tions about the Z-axis. We find the Hamiltonian that ro- 
tates the Bloch vector from its current position (as given 
by 9 and 0) through angles, 



56 c = N'(6t)- 
8<j) c = N\8t)- 



1ZEM 



tgo 

'izem) sin0 



(9) 
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Making a linear approximation and solving for the Hamil- 
tonian controls [H Xc and H Za ), we obtain expressions, 



H x 



h(se c 



2(8t) sn\(j) 

We) 



2(6t) 



(10) 



(In practice, although the expression for H Xc includes a 
1 / sin <f> term, removing the dependence upon <j) does not 
affect the performance of the guidance to a large degree 
and it dramatically reduces the bandwidth required for 
the control signal and makes the control system more 
resilient to time delays). 

Applying these controls requires manipulating the bias 
fields, A V x and Q x . Clearly there are limits to the size of 
the controls that can be applied using these bias fields. 
The two state approximation for the charge qubit is only 
valid as long as changes in the bias voltage are small 
|A14| <C 2e/C, so we impose constraints, such that 
|AT4J < 0.1 x 2e/C where AV Xo is the nominal voltage 
bias point and the voltage bias control fluctuates about 
this point AV X — AV Xo + AV Xc . Here the capacitance 
that is used to apply the gate voltage, which in turn 
controls the voltage applied across the qubit, is assumed 
to be the same as the qubit island capacitance. This 
is not necessary, changing the gate capacitance simply 
rescales the voltage bias and the associated behaviour 
under the action of the guidance/control. As long as 



this capacitance is known, the appropriate controls can 
be applied. For the magnetic flux bias, we must ensure 
that any fluctuation around the nominal bias point $ XQ 
is small enough so that the response is approximately 
linear. So $ x = $ Xo + A$ Xa where $ Xo = 0.25$ and 
|A$ Xc | < 0.05$o- Within this region, the cosine tun- 
nelling term is approximately linear in A^ Xr and we ob- 
tain the following relations for the control fields, 



A$ Xc = 



AV X . 



2H X 



nhv sin (^2- 



e 



(11) 



Beyond the limits given, the control is assumed to have 
saturated and any controls commanded are unachievable. 
This limits the number of states that are reachable from 
a given initial state, but as long as the desired state falls 
within the reachable set for the initial state, this should 
not be a major problem. For simplicity, we assume that 
the initial state corresponds to the ground state for the 
unperturbed Hamiltonian Ho at the nominal bias point 
(AV Xo and $ Xo — 0.25$o)- That is, for simplicity, we 
assume that the qubit has relaxed into the ground state 
via some (weak) dissipative process. Although it would 
be possible to prepare a different initial state by another 
process, such as some type of projective measurement 
interaction, this is not considered here. 

The guidance algorithm operates by integrating the 
evolution of the qubit for a small timestep and applying 
controls via AV X and <fr x that are determined by finding 
the ZEM angles from equations (6) and (7), converting 
these angles into an effective Hamiltonian using equa- 
tions (9) and (10), and finally converting these Hamil- 
tonian controls to bias values via equation (11). Apply- 
ing this procedure iteratively generates a time-dependent 
bias signal that can then be applied to a qubit and should 
provide the desired state at the desired time (t go — 0) as 
long as it is within the reachable set. Of course, since 
the control is currently open-loop, the actual state of the 
system is unknown. The controls are generated from the 
knowledge of where the state should be, if it started in 
the ground state and the controls had been correctly ap- 
plied. A wide variety of states may be prepared in this 
way and the controls have several distinct advantages: 
they have a comparatively low bandwidth, they operate 
via the standard bias fields without an additional exter- 
nal drive, and (for states within the reachable set) the 
controls vanish as t go — > 0. This last point means that 
the Hamiltonian is only weakly perturbed at the time 
when the state is required. Although the algorithm is 
presented as an open-loop control technique, it provides 
a natural generalisation to feedback (closed-loop) control, 
which is discussed in a later section. 

The difference between the desired state and the final 
(possibly mixed) state can be quantified in a variety of 
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standard ways. For the purposes of this paper, we use 
three common measures to characterise the performance 
of the guidance algorithm. The first measure is the fi- 
delity of the state, introduced by Jozsa ^U, which is one 
measure for how close the final state is from the desired 
state. The fidelity F for two density matrices, pf and pd 
is given by [l9|. 



F = F(p hPd ) 



Tr 



PdPfVPd 



which varies between zero and one (one being that the 
final state matches the desired state exactly). The second 
measure that we use in this paper is the trace distance D 

m, 

D = D {PfiPd) = \ Tl iPf - Pd] 

which also runs from zero to one (zero being identical 
states) and measures the separation of the desired and 
final states. The final measure is the distance from the 
surface of the Bloch sphere, which is a measure of the 
purity (or, conversel y, t he mixedness) of the state, and 
can be written as [illl6|. 



p = 2Tv[p 2 } - 1 

with pure states giving p — 1 and maximally mixed states 
giving p — 0. The purity is most important for mixed 
states, which are generated by the inclusion of stochas- 
tic terms, such as the measurement and feedback model 
discussed in a later section. 



IDEAL OPEN-LOOP PERFORMANCE 

Starting the qubit in its ground state, we can pick a par- 
ticular target state and calculate the control signal re- 
quired to generate this state at a later time using the 
technique described above. We start by examining the 
performance of the guidance algorithm when applied to 
an ideal system (with no errors) and then examine the 
performance in the presence of possible experimental er- 
rors. 

Although any target state can be chosen, an obvious 
candidate is the qubit excited state (which is a function 
of the nominal bias fields). As we have said, we assume 
that the initial qubit state will be the ground state, the 
qubit having relaxed into this state prior to the control 
being applied. Figure 2(a) shows the fidelity and trace 
distance achieved for a range of AV Xg and for three dif- 
ferent time periods t max — 20, 50, 100 qubit cycles and 
N' = 8. The purity is trivially one, since there is no 
stochastic evolution and the qubit state is therefore pure. 
We see from the graphs that the excited state is easier 
to reach (high fidelity and low trace distance) from volt- 
age values close to zero, AV Xo ~ 0. The fidelity is one 



for a wide range of bias voltages, even when the time 
available is relatively short, t max = 20 cycles. The main 
reason that the excited state is easier to reach near zero 
bias is that the ground and excited states at zero bias are 
equal (symmetric and anti-symmetric) superpositions of 
the two basis states, and the main controls that are re- 
quired are rotations about the Z-axis (<7 Z ), which are eas- 
ier to achieve than rotations about the X-axis or Y-axis. 
As the bias voltage is increased, the ground and excited 
states shift toward the poles of the Bloch sphere (and 
the qubit natural oscillation frequency increases, reduc- 
ing the time available), and more controls are required 
from the magnetic bias field (which generates a x terms) . 
However, allowing more time for the guidance (by ex- 
tending t max ) increases the ability to generate an ex- 
cited state, and for t max = 100 cycles the excited state is 
within the reachable set for the whole of the range AV Xo 
shown in Figure 2(a). Selecting the excited state as a 
target state has one more practical advantage. Because 
the exited state is stationary by definition, the produc- 
tion of an excited state is quite robust. Slight variations 
in the control fields or the time to go only generate small 
deviations from the excited state. 

Figure 2(b) shows the performance of the guidance al- 
gorithm for a target state that is not an energy eigenstate 
for any values of the bias fields. We select an equal su- 
perposition [6 = 7r/2) but with a phase 4> = 7r/4. Here 
we see that the reachable set (in terms of AV Xa ) is very 
much reduced compared to the previous example. Even 
though the fidelity is near one for a large range of bias 
voltages, the trace distance is significantly greater than 
zero until t max — 100 cycles and AV Xo ~ 0.2 x 2e/C. 
This is a result of the fact that the target state is not an 
energy eigenstate and is therefore non-stationary for all 
bias values. The control is therefore much more sensitive 
to small variations. 

Figure 3 shows an example of the evolution of the qubit 
state under the proportional navigation guidance and the 
control fields that were applied. The evolution starts in 
the ground state and rapidly spirals around the Bloch 
sphere under the influence of the control, and gradually 
approaches the excited state, spiralling in gradually as 
the controls applied reduce in size. This is a good ex- 
ample of the benefit of this approach, where the con- 
trols subside to zero as the system approaches the de- 
sired state. This has distinct experimental advantages 
because the bias fields will be static immediately prior to 
tgo = 0, which means that the Hamiltonian is not vary- 
ing rapidly when the desired state is required. It is also 
noticeable that both control fields contain a dominant 
frequency component that matches the coherent oscilla- 
tion frequency of the qubit, indicating that a coherent 
drive at the transition frequency is an important part of 
the control fields. 
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FIG. 2: (Color online) Fidelity (F - solid lines) and trace 
distance (D - dashed lines) achieved using guidance algorithm 
as a function of the voltage bias value for different values of 
tmax'- 20 qubit cycles (blue), 50 qubit cycles (green), and 
100 qubit cycles (red); (a) target state is the excited state 
at nominal bias voltage, (b) target state is an arbitrary state 
with = tv/2 and <\> = 7r/4. 



IMPERFECT BIAS FIELDS 

In classical guidance systems the guidance-control feed- 
back loop must be robust enough so that small pertur- 
bations from noise or imperfections in the control sys- 
tem are damped out and the system achieves its objec- 
tive. The main imperfections in classical guidance tend 
to come from uncertainties in the physical parameters 
that define the transfer function between the accelera- 
tions commanded by the guidance law and the actual 
accelerations achieved by the controls. In the qubit guid- 
ance case, this is more difficult because of the problems 
already mentioned in measuring the error signal and the 
algorithm discussed so far is an open-loop control system. 
Instead, we require that the fidelity and trace distance of 
the final state be weakly sensitive to small variations in 
the experimental control parameters. For example, the 



FIG. 3: (Color online) (a) Example trajectory of state on the 
Bloch sphere under proportional guidance with the excited 
state as the target state and A14 = 0.1 x 2e/C; (b) Control 
fields AV Xc (green) and A& Xc (blue) as a function of the t go 
for example trajectory in Fig 3(a). 



static (nominal) bias point may only be known to a cer- 
tain accuracy. We should require that the final state 
is approximately correct even if the static bias field is 
slightly off or if a small amount of dynamical noise is 
present in the fields. This is also of concern when the 
qubit is non-ideal, either the gate capacitance is only 
known to a finite accuracy (so the scaling of the voltage 
biases is inaccurate) or the Hamiltonian includes non- 
ideal terms (possibly due to variations in the junctions 
coupling the island to the bulk material). In these cases 
the controls applied will not necessarily generate rota- 
tions about axes exactly aligned to the the X- and Z-axes. 

In practice, we find that - in common with the per- 
formance of the guidance algorithm itself - the sensi- 
tivity to noise is dependent on the target state. For 
static errors of the order of A(AV^ Q ) ~ 10~ 5 x 2e/C 
and A(Q Xo ) ~ 10 _4< l>o, the performance of the algorithm 
for the example given in Fig 2(a) (i.e. the target state 
is the excited state) is very good. (The performance is 



7 



also good for time-dependent errors, as long as the cu- 
mulative errors during the control cycle are of the same 
order as these tolerances). Even for variations an order 
of magnitude larger than this, the performance is still ac- 
ceptable for the excited state over a comparatively wide 
range of bias voltages, even if the reachable set is signif- 
icantly reduced. The accuracy of the final state for the 
target state given in Figure 2(b) would be significantly 
less with errors of this size. The trace distance between 
the final state and desired state is noticeably larger when 
bias errors are introduced. For dynamical noise, we have 
also characterised the performance of the algorithm in the 
presence of white noise (i.e. uncorrelated with a uniform 
frequency distribution), and the performance is similar 
to that for static bias errors as long as the cumulative 
drift of the bias fields is less than the limits given above. 



FIRST-ORDER TIME DELAYS 

In addition to noise, an experimental system is also likely 
to contain other imperfections. The main one considered 
here is a constraint on the bandwidth allowed for the 
control signal. We use a simple method of restricting the 
bandwidth in the control system by introducing a first 
order time-delay, which acts as a low-pass filter and has 
a transfer function (Laplace transform) given by [2l|, 



where is the time delay constant. An example of this 
type of delay is a low-pass RC filter, which should be 
familiar from standard electrical circuit analysis. The 
effect of this filter is to exponentially damp rapid varia- 
tions in the controls, and in Figure 4 we show the effect 
of such time delays on the performance of the guidance 
algorithm. Time delays of the order of Td ~ 0.1 — 0.2 
cycles have little effect on the fidelity or trace distance of 
the final state for bias voltages close to zero, but toward 
the far right (higher voltage biases), the excited state be- 
comes harder to reach, indicating that the bandwidth of 
the control signal is larger at these extreme values. By 
contrast, the effect of first-order time delays on the case 
shown in Figure 2(b) would actually be less pronounced 
than that shown in Figure 4 because the range of bias 
values from which the target state is reachable is already 
comparatively small. 



CLOSED-LOOP PERFORMANCE 

In this section we consider the use of a simple mea- 
surement and Markovian feedback mechanism to demon- 
strate how the open-loop guidance approach could be 
adapted using existing quantum feedback techniques. As 
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FIG. 4: (Color online) Fidelity (F - solid lines) and trace 
distance (D - dashed lines) achieved using guidance algorithm 
as a function of the voltage bias value for different time delays: 
zero time delay (blue), 0.125 qubit cycle delay (green), and 
0.25 qubit cycle delay (red); target state is the excited state 
at nominal bias voltage and t max — 50 cycles. 

an example, we choose a simple model for photon detec- 
tion and (instantaneous) feedback. The model assumes 
that the qubit is weakly coupled to a lossy reservoir and 
that projective measurements are made on this reservoir. 
The results of the measurement are then used to modify 
the controls applied to the qubit. This model may not 
be entirely realistic, because of problems with detecting 
single microwave photons and with the large bandwidths 
required for a rapid feedback, but it demonstrates the 
general approach. The basic idea is to apply the guidance 
algorithm as described above and to modify the control 
pulses, to allow for the reduced time to go, when a photon 
is emitted and detected (detection is assumed to occur 
with efficiency rf). Where photons are not detected, the 
qubit will still be coupled to the lossy reservoir which 
will cause an additional (non-unitary) perturbation on 
the otherwise coherent dynamics of the system and this 
is allowed for in the modelling but not in the guidance- 
control algorithm. 

The measurement mechanism is modelled using a 
quantum trajectory approach |22L , corresponding to 
an unravelling of the Markovian Master equation for the 
qubit reduced density operator (after performing a par- 
tial trace over the lossy reservoir). In this paper, we 
choose the quantum jumps approach |22l l23l |. which is 
suitable for modelling spontaneous emission processes 
and is computationally efficient [24| • Physically, this un- 
ravelling corresponds to the detection or absorption of 
spontaneously emitted photons on a time scale that is 
significantly faster than any of the time scales present 
in the quantum system. All unravellings reproduce the 
Master equation evolution when averaged, and the indi- 
vidual quantum 'trajectories' for the qubit are described 
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by a model given in reference |24j . 

The spontaneous emission process and subsequent de- 
tection of the photon introduces quantum jumps that 
project into the instantaneous ground state of the qubit. 
The probability of a spontaneous decay occurring during 
a small - but finite - time interval, St, is 7 (c^c) St, where 

and c are the raising and lowering operators for the 
(instantaneous) qubit energy states respectively. During 
each time interval where no spontaneous decay occurs, a 
non-unitary evolution operator 

Cla(dt) = 1 - ~H qu St - ^cSt 

is applied to the qubit state. When a decay occurs, an 
operator 



Clx(dt) = 




is applied to the qubit state. Each run of the simulation 
produces a subjective 'trajectory'. For each trajectory, 
feedback control is invoked if a photon is detected, with 
probability rj, so that the state evolution is conditional 
upon the detections and then averaged over multiple re- 
alisations. The result is averaged over many runs to pro- 
vide an estimate of the mixed state density matrix for 
the qubit, from which we can calculate purity, fidelity 
and trace distance. (Some noise is still present in the av- 
erage density matrix and the performance measures, but 
this is relatively small and is due to the limited number 
of runs (typically 500-1000) which is dictated by compu- 
tational constraints). 

Figure 5 shows an example of the fidelity, trace dis- 
tance and purity for this simple closed-loop control sys- 
tem, corresponding to one of the examples shown in Fig- 
ure 2(a) with the damping rate 7 = 0.05/qubit cycle. If 
the damping rate were significantly smaller than this, the 
probability of a qubit undergoing a spontaneous decay 
during the period of control would be negligible and the 
system would reduce to the open-loop case already dis- 
cussed. If the damping rate were significantly higher than 
this, guidance-control (open-loop or closed-loop) would 
be impossible. For the case shown, the probability of 
undergoing a transition during the control period is sig- 
nificant. Without feedback (or very inefficient photon 
detection rj — - blue lines) the fidelity and trace dis- 
tance are very low for small voltage bias values, and the 
purity is fairly high. This is an indication that most of 
the qubits will spontaneously emit photons and decay 
back to the ground state. The main controls are applied 
near the start of the control period (see Figure 3(b)), 
which will tend to leave the qubit near the ground state 
once a photon has been emitted. In fact, the two prop- 
erties are related, since the qubit is only likely to emit 
a photon once the control has brought the system close 
to the excited state in the first place. As the bias volt- 
ages are increased, the excited state is harder to reach 
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FIG. 5: (Color online) (a) Fidelity (F - solid lines) and 
trace distance (D - dashed lines), and (b) Purity (p - cir- 
cles) achieved using guidance algorithm as a function of the 
voltage bias value for 7 = 0.05/cycle and rj = 0.0 (no feedback 
- blue), rj = 0.5 (green), and rj = 1.0 (red); target state is the 
excited state at nominal bias voltage and t ma x = 50 cycles. 



(and therefore occurs later in the control cycle) and the 
emission probability consequently goes down. This is the 
cause of the minimum in the purity for the rj = case, 
where the two effects balance out, so that the mid-range 
voltages are more likely to be mixed between the ground 
and excited states. 

Where the detection probability is non-zero, the green 
and red curves in Figure 5, feedback is allowed and the 
controls can be modified when a photon is detected. In 
the case of rj = 0.5 the purity is significantly reduced be- 
cause there is a chance of being near either state for most 
values of the bias voltage, either through decay and detec- 
tion, decay and non-detection, and non-decay. However, 
as desired, the fidelity increases as the detection proba- 
bility increases, indicating that the feedback is working 
correctly. Even so, even with 77 = 1, the closed loop per- 
formance does not reach the open-loop, non-dissipative 
performance. This is because the spontaneous emission 
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reduces the effective time available for the control, and 
multiple jumps are likely to occur for some bias voltages. 
The average number of jumps is dependent on rj and 
AV X0 , but for rj = 1 and AV X0 ~ two or more jumps 
are not uncommon. 

DISCUSSION 

In this paper we have presented a generalisation of a clas- 
sical guidance law to the problem of control of a qubit 
state on the Bloch sphere. We have chosen the propor- 
tional navigation guidance law because of its relative sim- 
plicity and its resultant widespread use in classical guid- 
ance and control systems. We have demonstrated that 
this guidance law can be used to generate an arbitrary 
quantum state from the ground state of a superconduct- 
ing charge qubit using the standard control fields (volt- 
age and magnetic flux bias). The controls produced by 
this guidance law are relatively robust to imperfections in 
the control fields and to first order time delays, implying 
that the control signals have a comparatively low band- 
width. This should make it possible to control the state 
of the qubit using signals fed through low-pass transmis- 
sion lines. We have also suggested a simple method to 
allow the generalised guidance law to be included in a 
closed-loop (Markovian) control system. 

Although the ability to control the state of a single 
qubit with a high degree of accuracy is important for pos- 
sible quantum information processing device, the ability 
to control the collective behaviour of multiple qubits is 
also of great interest. The ability to visualise the con- 
trol of the qubit state on the Bloch sphere is useful in 
understanding the guidance mechanism, but it is not es- 
sential. Generalising the guidance algorithm to higher 
dimensional settings (multiple qubits or N-level systems) 
simply requires an understanding of the group structure 
of the space and the ability to create appropriate control 
Hamiltonians from the generators of the group J25|. As 
long as the group generators (or the restricted set of gen- 
erators available to the control system) allow the state 
space to be explored fully, then it should be possible to 
generalise the guidance algorithm described in this paper 
to higher dimensional systems. 
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